A Simple Algorithm to Calculate S(n) 
by John C. McCarthy 


Introduction 


This short paper first outlines an "obvious" algorithm for 
calculating S(n) (the smallest integer m such that m! is divisible by 
n). Doubtless, there exist more subtle and efficient algorithms. I 


hope some readers will devise these and enlighten me concerning them 
through this journal. 


This is followed by a small scale investigation of the efficiency 
of the algorithm. 


Then there is a short discussion of a simple way of reducing the 
space required for storage of all S(n) for ranges of n. The storage 
space required for S(n) for all n which my routines can handle is 
considered. 


Heavily commented listings of an implementation of the algorithm 
in "C", sample output and timing data are included to help illustrate 
the algorithm. 


he Algorith 
The algorithm is described in detail at the start of the header 
file "“S(n).H”. Together with “S(n).c", this forms all the code 


necessary to implement the algorithm. Note that, for the S(n) function 
to work correctly, the function make_primes() must first be called from 
the main program. 


The code for printing S(n) and timing the routines has been 
omitted. These activities are both implementation specific and easily 
done. They are therefore left as an exercise for the interested 
reader. 


The algorithm hinges on finding the prime factors of n. 
Improvements on how this is done will most benefit its efficiency. 


To be practical, the given implementation of the algorithm only 
works for 0<n<2?3?. However, the algorithm is generally applicable to 
any non-null integer. 


Tables of S(n), constructed using the routines of "S(n).C", for 
the largest 2000 permitted n are included.’ My paging 
routines are rather elaborate. Using them (without printing!), it took 
2.4 hours to discover that 3,745,708 pages, as tightly packed as those 
shown, would be required to print S(n) for all O<n<23?. 


Efficiency of the Algorithm 


In a letter to R. Muller (about computing the Smarandache 
Function, July 19, 1993), Ian Parberry (editor of <SIGACT News>, 
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For the smallest 4800 numbers, see Ibstedt’s table (pp. 43-50) of 19 
this current journal. 


Denton, Texas) expressed that one can immediately find an algorithm 
that computes S(n) in O(nlogn/loglogn) time ('A Brief History of the 
"Smarandache Function"' by Dr. Constantin Dumitrescu, Department of 
Mathematics, University of Craiova, Romania). Disappointingly, a 
little analysis of the accompanying timing data on my TI85 advanced 
scientific calculator reveals that my algorithm is somewhat worse than 


this. 


Trying to fit the version 2 timing data to various O(f(n)), I 
obtained the following results (x=3355443200 and 10(0(x+99)-O(x-100)) 
is calculated for comparison with the last entry of the version 1 
timing data): 


10(0(x+99)-O(x-100) 


O(f£(n)) | Correlation | O(232-1) 
(milliseconds) 


Coefficient (years) 






















O(n) 0.9928879 0.6092 8909 
O(nlogn/loglogn)| 0.9944006 0.7906 11827 
O(n/n) 0.9997756 24.2 469178 


O(nvn) fits the version 2 timing data best, although the time it 
predicts for the last entry of the version 1 timing data is almost 3 
times too large. Hence, I assume the time complexity of my algorithm 
is a little better than O(n/n). 


As a rough upper limit on the time my program (on my 20MHz2 368DxX 
PC) would take to calculate S(n) for all O<n<2°2, let us assume that 
every value of n requires as much time as each n in the range of the 
last entry of the version 1 timing data (= 159111/199/10 = 79.9553 ms). 
In this “worst case", it would take 10.882 years. O(nvyn) time 
complexity predicts more than twice this value, which is a measure of 
how pessimistic it is. 


I would welcome a more rigorous analysis of the time complexity of 
my algorithm as I presently lack the necessary expertise. 


Simple Compression of Stored S(n) 


Without compression, each S(n) would be stored as a 32-bit (= 4 
bytes) value. Hence 2°? bytes (= 16 Gigabytes) would be required to 
store S{n) for all O<n<232. 


This requirement can be reduced considerably if we use the high 
bit of each each byte of each value to indicate if it is the last byte 
of the value. If the bit is set it means that further byte(s) are 
required and if it is reset it means that the byte is the last byte of 
the current value. This means that only 7 bits of each byte are used 
to form the numerical part of the value. Assuming that, as with Intel 
format, the values are stored low-'byte' {actually 7 bits) first, here 


are some examples: 


i) 127 requires seven bits and so just one byte (with high bit 


reset to indicate no further bytes). 
ii) 16,000 requires 14 bits. So it is stored as two bytes. The 
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first is 0 (16,000 mod 128) + 128 (to set the high bit 
indicating there is more to come). The second is 125 (16,000 
div 128) (with high bit reset to indicate no further bytes). 
This reads simply as 0 (with more to follow) + 128*125 (no 
more to follow). 

iii) A number stored as the three bytes 57+128, 93+128 and 125+0 
would similarly represent: 
57 + 93*128 + 125*128*128 = 2,059,961. 


The largest numbers that can be represented by a given number of 
bytes is thus as follows: 


127. 


I byte can code up to 27-1 = 

2 bytes can code up to 2?4-1 = 16,383. 

3 bytes can code up to 2?7'-] = 2,097,151. 

4 bytes can code up to 278-1 = 268,435,455. 

5 bytes can code up to 235-1 = 34,359,738,355 (or 8 times the 


largest unsigned long). 


For small values of n, the savings are considerable (400%). 
However, even large n often have small S(n). 


Using this technique to compress all S(n) calculated for some 
ranges of n (each range was also stored), I obtained the following 
results: 


range of n compression time taken size size after pkzip 
(seconds) 

1 without 4.5 40,008 19,836 

2,147,478,648 without 827.3 40,008 33,729 

-2,147,488,647 with 842.4 33,541 30,836 

4,294,957,296 without 1,066.2 40,008 34,320 

-4,294,967,295 with 1,085.1 34,330 31,634 


The results indicate that this compression is a little better than 
pkzip's (a commercial file compression utility). Application of pkzip 
to a pre-compressed file also gives a slight improvement. 


Assuming that the savings shown for the middle range of 10,000 n 
are the average of all ranges of 10,000 n, using my compression 
together with that of pkzip would permit storage of S(n) for all 
O<n<27?7 in about 3.0836*2*? = 12.3344 Gigabytes. So look out for sets 
of 19 CD-ROMs with all your favourite numbers on them! 


21st November 1993 
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S(n).H 


(c).1993.11.13.John.C.McCarthy 
"Sinin" 


Example Implementation of A Simple Algorithm to Calculate S(n), 
The Smarandache Function: 


Because there are more people familiar with C than with C++, this module 
has been written entirely in C (apart from "//" style comments). The 
module was compiled using Borland C++ version 3.1. 


For efficiency, n is constrained to the limits of an unsigned long. Hence, 
0 <= n <= 2°32 - 1 (= 4,294,967,295). (""" represents exponentiation). 
Although catering for n of vast magnitude is possible, it imposes heavy 
storage and processing overheads. The range of an unsigned long therefore 
seems a reasonable compromise. 


The algorithm depends on the most elementary properties of S(n): 
1) Calculate the STANDARD FORM (SF) of n: 
In SF: n = +/-(pl°al)*(p27a2)*...*(pr°ar) where pl, p2,...pr denote the 


distinct prime factors of n and al, a2,...ar are their respective 
multiplicities. 
2) S(n) = maxf S(pl7al),...,S(pr°ar)]. 
3) S(p°a), where p is prime, is given by: 
3.1) a<=p == S(p°a) = p*a. 
3.2) a>p == S(p7~a) = x < p*a. In this case, fortunately rare, x 


is the smallest integer such that p appears as a factor in the 

list of all integers > 1 and <= x at least a times. Let the no. 

of times p appears as a factor in the list of all integers > 1 

and <= y be fly, p). Then: 

f(y, p) = Ulintly/(p7i))] for i>0 while yo=(p7i). 

Hence, x is the smallest integer such that f(x, p)>=a. 

Note that between succesive integer multiples of p there are no 

integers which have p as a factor. The trick here is to look for 

the largest multiple of p (call it c), such that f(p*c, p)<=a 

(so that x = p*c, if f(p*c, p) = a, else x = p*(ctl)): 

3.2.1) c = a-2 (largest possibilty for c since f(p*(a-1), p)>=a 
when a>p (Note: f(p*(a-1), p)=a is not sought for slight 
performance gain)). 


32:2) 2 St Coe x pla 

3.2.3) While(2>a): 
3.2.3.1) d = no. of times p appears as a factor of p%*c 

= (no. of times p appears as a factor of c) + 1. 

3.2.3.2) c = c-l (next largest possibility for c). 
3224563) 2 = 2-4 (= f(p*c; DJ- 

3.2.4) If(z2<a), x = p*(ct#l). 

3.2.5) Else X = p*c. 


To calculate the prime factors of all 32-bit n requires use only of primes 
< (2716) (i.e. all primes expressible as an unsigned short integer). This 
is because any factor of n remaining after division of n by all its prime 
factors < (2716) is simply a prime. Since there are only 6542 16-bit 
primes, the program first creates a list of these (which only takes about 
4 seconds on my 20 MHz 386DX PC) so that they never have to be 
recalculated, thus saving much time. 


S(n).H 


#define PRIMES16 6542 // The number of 16-bit primes 

#define MAX_SFK 9 /* max. distinct primes in the SF of n. The smallest 
number with more than 9 distinct primes is the product of the 10 smallest 
primes (= 6,469,693,230), which is substantially more than the largest 
integer expressible as an unsigned long. Hence, 9 distinct primes are 
more than ample. 

* 

typedef unsigned long u_long; 

typedef unsigned int u_int; 

typedef enum {false, true} boolean; 


struct SF struct { 


int sik; // no. of distinct primes 
u_long sfp[(MAX_SFK]; // the distinct primes 
int sfalMAX_SFK]; // respective multiplicities 


b; 


extern u_int prime[PRIMES16+1]; // list of all 16-bit primes 
// plus terminating zero. 


void make_primes(void); // construct list of all 16-bit primes (prime[]). 
= // Must be called before calls to getSF() or S(). 

void getSF(u_long n, struct SF_struct *SF); // calc. SF of n and store in SF 

u_long S(u_long n); // calc. S(n) 

u_long Spa(u_long p, int a); // calc. S(p°”a) where p is prime 

int f(int x, int p); /* the number of times the prime p appears as a factor 
in the integers from 1 to x inclusive. This function is only called from 
Spa(p, a) when a>p with x=p*(a-2) (refer to item (3) of algorithm outline 
above). Max value of (a) occurs when p is a minimum, n is a maximum and 
(p°a)=n. So, (27max(a))=max(n)=(2732)-1. Hence max(a)<32. So, x<60 
when (a) is at its max. Max value of p (and x) occurs when a=ptl and 
(p°a)=max(n). So, max(p)7 (max(p)+1)=(2732)-1. The upshot is that 
max(p)=9 when a=10. Hence, max(x)=72. This explains why it is safe for 
x, p and the return value of f(x,p) to be passed as ints. 

ai 
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S(n).c 


/* (c).1993.11.13.Jonn.C.McCarthy 


= / 


S(n).c” 


Example Implementation of A Simple Algorithm to Calculate S(n), 
The Smarandache Function: 


This is the code for the module. Refer to "S(n).h" for details. 


#include “S(n).h" 


u_int prime(PRIMES16+1]; // allocate storage for list of all 16-bit primes 


// plus terminating zero. 


void make_primes(void) 


{ 


} 
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u_int *pp; // ptr to last prime so far of prime list 
u_int *tp; // ptr to current test prime 


u_int p; // number being tested for primality 

pp=prime; // point to start of prime list 

*pp=2; // set first prime to 2 

*++pp=3; // set second prime to 3 

p=5; // next possible prime. N.B. p is kept odd so that trial 


// division by 2 is unnecessary. 


while(true) { // infinite loop!: 


} 


tp=prime+l; // point to first odd test prime 
// whilst test prime <= Vp: 
While(((long) *tp)*(*tp)<=p) { 
if(!(pé*tp)) { // If current test prime divides (is factor of) p: 


Ppt=2; // try next odd number 
if(p<*pp) { // done when p overflows: 
*++pp=0; // terminate list 
return; 
tp=primet+l; _// point to first odd test prime 
else ++tp; // Else point to next test prime 
} 
// no prime <= Vp divides p so p must be prime: 
X++pDp=p; // so store it next in the list 
pt=2; // try next odd number 
if(p<*pp) { // done when p overflows: 
*++pp=0; // terminate list 
return; 


} 


S(n}).c 


void getSF(u_long n, struct SF_struct *SF) 


u_int *pp; // ptr to current prime 
u_long rT; // ‘residue' of n remaining for factoring 


SF->sfk=0; // no. of distinct prime factors discovered 


r=n;} 


pp=prime; // point to start of prime list 


// whilst current prime <= /r and prime list not exhausted: 


while(((long) *pp)*(*pp)<= 


if(!(ré*pp)) { 
SF->sfp(SF->sfk]=*pp; 
SF->sfa(SF->sfk]=1; 
r/=*pp; 
while(!(r%*pp)) í 
SF->sfalSF->sfk] ++; 
r/=*pp; 


SF->sfk++; 
} 
++DDj 
} 


if(rm1) { // 
SF->sfp(SF->sfkl=r; 
SF->sfa(lSF->sfkJ=1; 
SF->sfk++; 
} 
} 


r && *pp) { 
// if current prime is a factor of r: 
// store current prime as next prime of SF 
// set its multiplicity to i 
// ‘divide out' current prime 
// while current prime factors r: 
// increment multiplicity 
// ‘divide out' current prime 


// increment count of distinct prime factors 
// next prime 
If n contains prime > 2716: 
// store it as last prime of SF 


// set its multiplicity to 1 
// increment count of distinct prime factors 
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u 


{ 


} 


S(n).C 


long S(u_long n) 


struct SF_struct SF; // to store SF of n 


int sfi; // index of current term of SF of n 

u_long Sn; // current guess at S(n) 

u_long x; // S(current term of SF of n) where it might exceed 
// current value of Sn. 

if(n==1) return 0; // special case 

getSF(n, &SF); // cale. and store SF of n 


// First guess at S(n) is S(p°a), where p is the largest prime in the SF 

// of n and a is its multiplicity. This pre-empts the calculation of S(p^a) 
// for the remaining terms where, as is likely, p*a for these terms is <= 

// this initial guess (since S(p°a) <= p*a always): 

sfi=SF.sfk-1; 

Sn=Spa(SF.sfp{sfi],SF.sfa(lsfil]); 


while(sfi>0) { // while more term(s): 


sfi--; // next term 
if(SF.sfp[(sfil*sFr.sfalfsfiJ>Sn) { // if this term may have larger S(p’a): 
Xx=Spa(SF.sfp{sfil,SF.sfa(sfil); // calc. it 
if(x>Sn) Sn=x; // if new max., update Sn with it 


ry 
} 
return Sn; // That's all folks! 


u_long Spa(u_long p, int a) 
{ 


} 
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// Refer to item 3) of the algorithm description in S(n).-h. 

int c; // largest multiple of p such that f(p*c, p) <= a (eventually!) 
int z; // E(p*c; p) 

int m; // used to calc. no. of times p appears as factor of c 


if(a<=p) return p*a; 


c=aąa-2; 
z=f(p*c, p); 
while(2>a) { 
// d in items 3.2.3.1) and 3.2.3.3) of algorithm description is implicit 


// nere: 
— 
M=C--; 
while(!(m&p)) { // while p divides m: 
— 
m/=p; // ‘divide out' factor of p fromm 


} 
} 
if(z<a) return p*(c#+l); 
else return p*c; 


int f(int x, 
{ 
int k=0; 
int xdp; 


xdp=x/p; 


S(n).c 


int p) 


// count of appearance of prime p as a factor in the integers 


// from 1 to x. 
// successive divisions of x by p 


while(xdp>0) { 


Kk+=xdp; 
xdp/=p; 


return k: 
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SMARANDACHE FUNCTION, S(n). for n=#4294965296 to n=4294967235 


\S of A plus: 


AN 6 1 2 3 4 s 6 
4294965296 268435333 165837 2799847 60492469 14316551 3226871 49941457 
4294965306 $531 4294965307 1194373 101377 5174657 66569 967 
4294965316 1073741329 4583741 102261079 273564667 1602599 4919777 2147482681 
42949653286 2147482663 648191 7255009 26712 537427 4294965331 3741259 
4294965336 21517 11392481 $520521 2199163 12632251 40139863 2286989 
4294965346 3919 46294965347 6911 90497 83641 14631655117 1228537 
4294965358 14969 158397 3109 1431108907 2339 4294965361 2147482681 

49254 655123 $5137 1897411 35 7 
4234963376 5527 330587 154573 273857 RETETE 20452216! 7333103 
4294365386 2147482693 21367987 3567247 30034723 7535027 35449 268435337 
4294965396 357913783 4040419 35461 27012361 1952257 377513 70867 
4294965406 235289 2347 154807 529393 5167 14461163 672349 
4294965416 5179 625997 22138997 4013 23860919 15107 17747791 
4294965425 816223 9965117 12109 6040739 321239 28793 5772803 
4294985436 65983 1530091 79538397 2473 13421787 1551089 308857 
4294985446 178793 655793 36997 564697 28633103 11732 121313 
4294965456 3329 4294965457 587067 423191 70019 4294965461 415337 
4294965466 10275037 14364433 7615187 1822217 49699 88873 169681 
4294965476 157 1194041 21474827393 1935541 18233 22031: 718943 
4294965486 102281083 4294965487 14639 973253 18673763 2754949 7253 
4294965496 50951 14387 715827583 390451409 1227133 835759 16909313 
4294965506 17903 72091 29483 3304221 255197 2017 193583 
4294965516 357913793 99882919 4699087 12236369 53687069 832519 715827587 
4296965526 2147482763 252645931 57193 4294965529 384509 1903391 32789 
4294965536 5335551 28881 88919 3296213 2468371 6529 230664! 
4294985546 19979 107077 1073741387 23339 23216983 622549 89478443 
4294365556 4948117 20355287 1038937 12569 28057 2293 7561 
4294965586 13174741 1431655189 67108837 397057 20452217 13709 51829 
4294965576 16268809 16183909 3081037 1431655193 12302541 4294965581 10374813 
4294985586 19701677 139 2843 2549 429496559 97392151 8301159 
4294965596 1073741399 218941 961763 226050821 198841 4294965601 2147482802 
4294965606 154373 118361 53687070% 923053 3283 35527 307 
4294965616 515231 6294965617 2957 89611 $147 7495577 2147482811 
4294965628 682513 1828423 46684409 2459 251699 9137 170327 
4234965636 9177277 682933 2147482819 1431655213 2551 4294965641 784039 
4294965646 2937733 209623 83478451 1197467 85899313 82245873 100547 
4294965656 20731 26371 147827 423934965659 1745922 16582879 118757 
4294965666 9061109 17111417 119557 4603393 629496567 4294965671 59652301 
4294965676 1073741419 28825273 715827613 6294965679 191873 18013 505409 
4294965686 195225713 1431655229 400949 2038427 47721841 4294965691 130262 
4294965696 2729 58741 1753 477218411 42949657 3373 94399 
4294965706 1059439 594953 1151 75431 61356853 13267 72727 
4294965716 220889 1286303 2147482859 889412 1556147 4294965721 118273 
4294965726 1716613 8493912 601873 8783161 79099 494413 23057 
4294965736 536870717 4294965737 179 7517 16519099 130150477 1589601 
4294965746 2147482873 887573 63161261 650851 121843 226050829 29173 
4294965756 2818219 4294965757 7877 204522179 1677722 18755308 89669 
4294965766 10269031 4294965767 Z3it 24793 7229 477218419 1000691 
4294965776 268435361 14046961 1162687 603989 4933 12377423 836573 
4294965786 4286393 244157 84407 1811 3539 41638697 7433 
4294965796 1337 $5823 57649 919 21474829 7052489 102559 
4294965806 23831 41357 24317 1160860157 143165527 38008547 19627 
42949652816 $761 677119 349241 2165833 57527 4294965821 29783 
4294965826 49829 2129 357913819 4799 38669 651641 535870729 
4294965836 359231 7883 5355319 4294965839 1626861 4294965841 2147482921 
4294965846 1918247 4294965847 36599 1431655283 721843 9546439 32749 
4294965856 7129 237173 72211 148192271 216748293 11098103 309391 
4294965866 107069 1186127 1073741467 33353 407879 613566553 2638435367 
4294965876 63113 36749 306783277 477218631 107376147 186737647 28031 
42949658856 2147482943 4294965887 2179 2081903 2161 366997 6869 
42949658596 19289 84121 2147482949 613566557 24316553 252645053 2147482951 
42934965306 6871 364691 10673741477 110127331 429496591 4294965911 87767 
4294965916 12064511 75217 187759 1848: 103643 405683 25799 
4294965926 256049 3119 1574401 52181 161913 97553 831713 
4294965936 89478457 4294965937 1821 130150483 214748297 490603 839703 
4294965946 74051137 263 3933119 42949659349 4651 79259 3947579 
4294965956 5417 33693387 6221 138547289 627919 90034727 306783283 
42949653966 1192061 4294965967 2591 8699 991969 4294965971 13597 
4294365976 76695821 4294965977 683893 1322 690509 1431655327 911 
4294965986 7430737 36709112 26188817 15259 229 61333 117709 
4294965996 6907 5956957 33331 42389 131 4294966001 715827667 
4294966006 36398017 4294966007 178956917 60492479 13159 243765 1073741503 
4294966016 15101 14737 415133 2573 $939 40031 19701679 
42943966026 279511 35495587 258173 254879 2225371 116080163 1439 
4294966036 $6512711 2235797 65075243 47197429 1105812 329927 93368827 
4294966046 238001 163003 134217889 273617 9544369 3256229 3343 
4294966056 4519 2213 2247483028 $9069 2788939 54366659 24683713 
4294966066 2147483033 613566531 3224449 9695183 10475527 3853 949223 
4294966076 103613 477218453 182129 91382257 1033 12521767 303617 
8234966086 79536409 42949656087 $701 1431655363 4013987 11392483 357913841 
4294966096 20543 44278001 45137 4294966099 829 85733 306783293 
4294966106 2147433053 669937 1973741527 35083 233549 4554577 7064089 
4294965116 131441 3271109 2147483059 33294311 22171 4294966121 14035837 
4294366126 2749 73583 257863 4294966129 674249 1759 242873 
4294966136 536870767 $21 2147483069 252645067 1213 1328477 58040083 
4294966146 1993349 340519 83873897 9697 40009 8641783 4153 
4294966156 4358559 186737659 77647 390451463 53687077 5168431 288679 
4294966166 15675059 15580% 905347 251 13015049 82139 16067 
4294966176 481067 4294966177 74051141 2237 9336883 390451471 55063669 
4294966186 126613 4294966187 18837571 80737 119737 21367991 24493217 
4294966196 1012009 47279 2008871 7548271 1022611 352133 6317 
4294966206 1063637 2113 5182219 57847? 1633067 2245147 $25571 
4294966216 14510021 4294966217 43481 237409 10853 1748053 30246241 
4294966226 16097 890887 $787 446323 199673 4294966231 536870779 
4294966236 1121987 4294966237 95287 1199071 246271 34239 23743 
4294966246 2147483123 566693 462421 28676809 3307 1431655417 1109 
4294966256 383479313 477218473 171491 12744707 13463 5501 16393001 
4294966266 163543 3$5573 1073741567 381673 61356661 96263 683 
4294966276 859681 87652373 7215827713 2363 9781287 46073 6691 
4234966286 165193011 15391 143779 16417 2563 311 632161269 
4294966296 304867 4294966297 3883333 332711 26141 72796039 389557 
4294966306 251197 2458481 357913859 1276653 919693 477218479 2753 
4294966316 26188819 5940479 103191 6074917 5951 226050859 139891 
4294966328 1571 116080171 5368767391 159972827 430789 19489 174337 
4294966336 17159 4294966337 1411889 325697 1614649 2706343 2147483171 
4294966346 195225743 68174069 77813 39403361 28633109 135347 5209 
42943966356 56893 8761 2147483179 75350287 107053 146681 23091217 
4294966368 283571 4294966367 581029 401887 429496637 1431655457 117877 
4294966376 26249 1907 113025431 390451489 14251 908603 108827 
4294966386 3851 1867376639 151637 2464123 273043 33818633 $90333 
4294966396 32783 1013681 17459213 399941 1342177 $557 48761 
4294966406 216098 1431655469 6047 41698703 3497 22067 97612873 


7 
231547 
4294965313 
22967729 
73577 
335613 
2545919 
23719 
104755253 
4294965343 
3749619 
26349481 
4294965413 
294522163 
746561 
2617 

13997 
138547273 
4001 

84869 
30509 
837519 
70383 
142799 
6381821 
21831861 
179623 
77933 
20507 
12781 
72796027 
33294307 
42794965613 
330381971 
1302689 
7027 
6400843 
28071671 
4294965673 
4294965633 
45613 
13054607 
2238127 
249287 
4294965733 
7309 
142837 
82183 
19792469 
42223 
4294965793 
556559 
14312655271 
81037091 
1395829 
48661 
346117 
42364061 
65173 
79309 
620923 
2633 

7823 
125941 
2243 
1290167 
23718358 
1431655321 
311749 
74471 
53282861 
1322749 
471301 
63372395 
3637 
4294966043 
238649 
351263 
4294966073 
46182431 
1081583 
340573 
1672695 
ZE33333 
665783 
732233 
4294566153 
4294966165 
6661 
817933 
14965039 
18313 

6719 

1361 
1739557 
4294966243 
116086163 
23671 
2217329 
35537 
6721387 
33747 
390451483 
49367429 
613566619 
252645079 
54133 
1762399 
4294966373 
1390728293 
1057613 
199001 
149333 
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3 
337867 
113025403 
2117833 
858893 
379691 
715827559 
118423 
641 
6628033 
2147482697 
511549 
34301 
39727 
59003 
28871 
2465537 
336870683 
5549051 
1073741371 
28661 
294337 
306783251 
10737413812 
12132671 
64319 
217599 
219304599 
298801 
29498393 
2970239 
40459 
78487 
4409 
2147482817 
4549 
160033 
20929 
93368819 
1360889 
113025413 
291907 
1831 
153391633 
2247482867 
6701 
2147482877 
1073741441 
685075239 
536870723 
306783271 
1431 
17471 
20719 
238609213 
1757351 
4393573 
23565273 
2147482937 
1091 
6949783 
1999 
7388279 
352931 
43876183 
536370743 
1494421 
241453 
1390253 
$1867 
S276371 
33053 
2318283 
529981 
333691 
3579133837 
126322531 
28627 
11387 
7110871 
8893 
28957 
6927364? 
1823 
31122943 
4793489 
2147483077 
578213 
9629969 
97987 
191 
2362623 
21303 
$263439 
4029049 
1073741561 
763957 
821 
22147483137 
58169 
74051143 
1234371 
13230377 
1600211 
49941469 
178956931 
209331 
1367 
1790557 
268435339 
126322541 
357913867 
1679033 


9 
2798023 
286331902: 
24342659 
858393067 
4621 
1493 
38629 
1307 
858993077 
86357 
977239 
24683 
1961054801 
14353457 
858993089 
2207 
168331 
6079 
8363 
137393 
72679 
18841? 
1847297 
122713301 
106747 
891997 
78090283 
1931 
337 
59529007 
5335361 
21859 
3061 
3193283 
401453 
46961 
858993133 
302:1 
12347 
7880671 
286331047 
858993143 
29153 
859853 
1553333 
6270023 
4457 
65497 
16943 
286331053 
8753 
122713309 
19088737 
4999 
$9409 
733 
7933 
15959 
3843491 
78090289 
858993181 
344387 
171798637 
3886847 
19927 
122713313 
8589931393 
457 
80513 
$383993:1393 
32497 
1340083 
19219 
11603 
1123769 
277383 
16231 
171798643 
68791 
13634813 
1433 
78090233 
11433243 
2340581 
27709459 
286331077 
122713319 
171798647 
86531 
365063 
19976587 
28633108: 
15618059 
812671 
15070057 
9323 
29620457 
179 
856993257 
37347533 
40904441 
313547 
197243 
397 
8538993269 
858993271 
25733 
701219 
73727 
19991 
838993281 
858993283 


SMARANDACHE FUNCTION. 


Stn), 


\S of n plus: 


nA 8 
4294966416 1877 
4294966426 13564081 
4294966436 37571 
4294966446 391%627 
4294966456 26861 
4294966466 306783319 
4294966476 38891 
4294966486 63273653 
4294966496 7833159 
4294966506 234467 
4294966516 306872 
42394966526 22351 
42949366536 8521759 
4294966346 165191021 
42949665546 73883 
4294966566 34781 
4294966576 1405421 
4294966586 389049089 
4294966596 383 
4294966606 13661 
4294966616 $7013: 
4294966626 79536419 
4294966636 8738671 
4294966646 21474833253 
4294966656 17449 
42949366666 556487 
4294986676 757 
4294966636 17939 
4294966696 5336870837 
4294968766 2147483353 
4294966716 119304631 
4294966726 S$Z2377643 
4294965736 1920667 
42943966746 115811 
4294966756 1073741689 
4294966746 321047 
4294968776 1844917 
4294986786 1490273 
4294966796 1789 
4294966806 761 
4294966816. 3257 
4294966826 29417581 
4294966836 48623 


2532388938 7147317453 


42949866866 191449 
4294966878 1073741719 


4294966886 306783349 
42949663896 1489 
42949663996 3639 
42949663916 7307 
42949663926 715827822 
4294966936 18121 
42949669456 20069939 
4294966956 1381907 
4294966966 70067 
42949663576 67108859 
4294966966 12358383 
42949663936 119159 
4294967006 114838693 
4294967016 34149 
42934967026 23917 
4294967036 2445881 
4294967646 16111 
429496706056 $347 
4294967066 15761 
4294967076 119304641 


429493967086 2147483543 


4294967096 76695841 
4294967106 1583 
42949671156 83841 


4294967125 2147483563 


4294967136 13147 
4294967146 933688351 
4294967156 1673741789 
4294967166 1262483 
4294967176 87481 
4294967186 152293 
4294867136 2521 
4294967206 1809168 
4294967216 37807681 
4294967226 22637 
4294967236 417961 
4294967246 28111 
4294967256 59652323 
4294967266 367531 
4294967276 361651 
42949672836 715827681 


1 
813566631 
4294966427 
17674739 
4294966447 
288427 
130150499 
4294966477 
613566641 
49367431 
416623 
49891 
477218503 
482993 
252545091 
1129957 
$3563 

233 
75350291: 
444i 
538835159 
105373 
45751 
3390382949 
12816397 
42934966657 
4234966667 
8251772 
390451517 
16582883 
975907 
3607 
64103981 
30460757 
13041 
58118863 
191 

9133 
148102393 
43383503 
99882949 
1516049 
13229 
613566692 


50734883 


999 
4294986877 
338213 
330382069 
1782607 
169286698 
4294966927 
2403451 
5990191 
3536998 
226050893 
163487 
569851 
4294966997 
65519 
795217 
18541 
1413283 
39153 
252645121 
132987 
23761 
4294967087 
1431655699 
2135737 
10058471 
21841 
138547327 
13990121 
1932059 
3229 
1047535297 
204522247 
4294967137 
21382281 
1431635739 
252645131 
825473 
65543 
$2613 
37549 
1861711 
65339 


for n=4294965296 


2 

48677 
7017919 
809 
1172207 
6567227 
2848121 
3308911 
1789586937 
2147483249 
1153 
1951 
26777213 
2147883269 
16193 
3919719 
463219 
37633 
1331693 
97171193 
29828157 
764301 
12064523 
7839 
1980223 
4121849 
152239 
49941473 
2383 
21691751 
4492643 
9193 
13765919 
11864349 
33653 
1299143 
27751 
99223 
34019 
2147483399 
32173 
3539171 
985302139 
15449521 


284233 


1073741717 
5976793 
23633 
1246363 
357913909 
82339 
372357 
42107519 
58512723 
23087 
59852319 
2147483489 
602887 
26107 
5835553 
14412641 
32537629 
8355967 
141319 
83459 
2599859 
52390571 
1132639 


858639 
1467863 
2251031 

2147483579 
180: 
317159 
3449 
12966907 
18539 
52377649 
10631107 

29389 

1877172 

2147483629 
29587 
18946081 
2089 


to nz=42$4967295 


3 

8088449 
$5369 
15913 
1137137 
172787 
17970571 
7001 
3455323 
4999961 
477218501 
124133 
613586647 
16086017 
39023 
339382943 
$11123 
138547309 
12739 
688627 
22233713 
42943665619 
23469763 
4294966839 
284831 
139381 
214309 
18979 
36709117 
209961 
16976153 
20164163 
208889 
32232983 
14316355383 
198391 
4234966769 
877218531 
85199 
17821438 
294522229 
1766749 
4294966829 
83509 


3431333 


2909 
911689 
1931 
154423 
4294966908 
641327 
8263 
186737693 
47197439 
397351 
97189 
250791 
1588963 
16339673 
1597 
1569053 
4294967028 
295493 
477218561 
1755197 
2408843 
110127361 
613566727 
96553 
1431655703 
15053 
161471 
159072857 
2944249 
5156023 
179383 
62929 
4294967189 
46182443 
48869009 
113593 
69779 
70409299 
15966421 
4817 

35809 
4294967279 
21039 


4 

251 

2939 
261251 
965161! 
2029021 
143165549 
681 
429496643 
86767 
81083 
8259551 
443237 
396947 
12271333 
64373 
201547 
214748329 
47721851 
32587 
1523 

6551 
1297573 
5821 
2202547 
86627 
42949668687 
63953 
Z21401 
37119 
2163165557 
431: 
144383 
71382779 
17179867 
4738 
47721853 
16519103 
7459 
3579139 
429496681 
214748341 
35081 
349753 


23066334 


233 
28843543 
3491843 
6135667 
2029 
2154562 
53436667 
19522377 
33923 
36497 
7621 

347 

98081 
4294967 
143165567 
180007 
11608019 
24% 

3009 
214748353 
1431655869 
1607374177 
429496709 
362 
1252177 
236507 
463319 
214748357 
85899343 
647 
184571 
34981 
143165573 
31033 
25264513 
25411 
13834733 
23327 
4111 
32771 
113861 
127 
22605091 


5 

30899039 
204522211 
4294966441 
2669339 
19611719 
12763 
330352037 
159072833 
813566643 
1864189 
345727 
21210873 
2207341 
1431635517 
119083 
23577 
21377 
42943966591 
465781 
119127349 
1693 
312839 
15671 
4294966651 
4294966661 
1237 
24826397 
109451 
75350293 
2280917 
2194 

84349 
40129 
4141723 
$9417 
116089183 
2387419 
362979 
93229 
1072937 
13899569 
91382273 
38331 


902878 


330332067 
143186855627 
8312 
113783 
1431655837 
36092159 
104755291 
8297 

9817 
447439 
1431655657 
4294966981 
3427747 
116127358 
22453 
365747 
S683 
10716641 
21582749 
7121 

85551 
417839 
21467 
181967 
4294967111 
102871 
128399 

1163 

77249 
4294967161 
6518339 
3383 

82791 
11576731 
477218579 
39403369 
4294967231 
1431655747 
1664071 
281 

3449 

8951 
4294967291 


8 

4804213 
377 
614093 
78301 
126322543 
8468323 
321143 
$6512717 
5413 

443 
3984199 
4493 
715827737 
5336870813 
52377841 
31907 
22139003 
222353 
293693 
5088823 
2126221: 
719 
74051149 
1073741663 
20063 
8659297 
86477 
357913891 
38569 
336870839 
482689 
-397 
195225761 
14913679 
46411 
107021 
715827797 
536870849 
93503 
32779 
26561 
3499281 
29347 


1030393 
1345541 


389040033 
1073741723 
725827817 
790051 
8555711 
37943 
306783353 
$36870869 
217643 
449 

202397 
323027 
4106087 
21913097 
5549053 
536870879 
96703 
357913921 
4097 
1226181212 
192261121 
3209 
66377 
182423 
12853183 
LO73 741783 
2476913 
16381 
34919 
357913931 
195225781 
23342213 
2385609289 
292493 
813748 
8191 
145681 
1527371 
19867 
5368709909 
733659 
331 


7 
39045149393 
389849 
2431655481 
3929521 
99882941 
1949 
21768712 
6806603 
1431655501 
390487 
9565627 
32611 
613566649 
4294966553 
617359 
1862:63 
4294966563 
21367893 
208261 
7993 

$2347 

8597 

14349 
799361 
756023 
1033189 
12030719 
54366667 
138547313 
1431655571 
63371 
889043 
5333459 
69463 
15176561 
24265343 
6338557 
3390382061 
1431655601 
4294966813 
299078939 
477218537 
192599412: 
rs783 
35983 

1327 
138653 

911 

59791 
13687 
453199 
127447993 
4294966943 
1431655651 
87652387 
390451543 
1431655661 
18341 
554977 
375467 
17402 
3764213 
9851 

372012: 
768193 
16086019 
390451353 
11332367 
68174081 
88463 
17322 
80021 
4294967143 
33813639 
4886197 
$733 

14767 
29333 
26029 
735167 

3089 

3357 
619159 
1431655751 
64103989 
395959 
1226783 
464773 
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8 

7213 

12541 
41539 

199 
24855313 
2147483237 
1973741621 
7866239 
336870813 
1256573 
10631 
29417579 
149713 
238609253 
330767 
$48199 
6170923 
298303 
1073741651 
9061111 
11719 
306783331 
153151 
195225737 
63997 
715827779 
1073741671 
1028981 
12782639 
701 
29239277 
4079 

362909 
79841 
756689 
306783341 
178481 
10683997 
4481 

4021 

25679 
27059 
275389 
sIn 
47731 
611993 
2753148439 
50647 
22133 

2797 

3271 
335204847 
9151 
2147483477 
1673741741 
715827829 
78341 
2147483497 
1325690712 
12707903 
141281481 
3392543 
34636832 


342283 
21692753 
$197 
675097 
81031 
186689 
7321 
70853 
3817907 
2147483587 
6197 
7134437 
97612891 
653981 
1726273 
318949 
45083 
862097 
262657 
6118187 
46684427 
2147483647 


9 
37266219 
636763 
65622! 
5039 
858993293 
1338: 
2951867 
858993299 
6151 
26407 
7469507 
40123 
56509 
32443 
16207421 
364753 
9439487 
179593 
286331107 
490573 
16901 
2219621 
4567 
15877 
26030101 
120223 
5867 
288331113 
218629 
2111 

4241 
924643 
838993349 
33671 
277069463 
15618061 
134523 
78311 
317323 
88129 
612951 
8589933967 
1579 


65557 
4517 


22930649 
221219 
3651611 
95443709 
858993383 
36291 
218398 
858993389 
683 
393855 
12113 
41177 

193 
§458597 
358993403 
10141 
858993407 
11767033 
1228889 
1997659¢t 
791699 
12917 
395303 
91983 
2863312141 
14753 
129347 
$76119 
515293 
838993433 
7283 
44171 
435210181 
286331147 
13109 
231223 
22025473 
69001 
4259 
95443717 
171798691 
122713351 
65537 


Timing of S(n).c Module for Calculation of Smarandache Function, version 1 





Time taken to calculate S(n) depends on how easy it is to factor n. Less 
time is required if n has "small" prime factors. So, in the following table, 
the values of n shown are the mid-points of ranges (n-99 thru n+99). Times 
shown are for calculating S(n) for all integers in each range 10 times over: 


n time (ms) 





100 268 

200 308 

400 345 

800 387 

1600 432 

3200 490 

6400 571 
12800 661 
25600 766 
51200 919 
102400 2450 
204800 4036 
409600 5670 
819200 7977 
1638400 10423 
3276800 13004 
6553600 16302 
13107200 23438 
26214400 29642 
52428800 37011 
104857600 50330 
209715200 62363 
419430400 77888 
838860800 108179 
1677721600 158480 
3355443200 159111 


30 


Timing of S(n).c Module for Calculation of Smarandache Function, version 2 


"Time to n” is the time taken to calculate S(n) for all n <= that shown. 
"Time add.” is the time taken to calculate S(n) for all n> previous n and 
<= current n. All times are in milliseconds (as per version 1): 


n Time ton Time add. 


50000 18223 18223 
100000 66763 48540 
150000 139191 72428 
200000 229634 90443 
250000 335252 105618 
300000 452539 117287 
350000 579419 126880 
400000 715146 135727 
450000 859963 144816 
500000 1012335 152372 
550000 1171221 158886 
600000 1336899 165678 
650000 1508825 171927 
700000 1686808 177983 
750000 1870023 183215 
800000 2058983 188961 
850000 2252457 193473 
990000 2450892 198435 
950000 2653620 202728 
1000000 2860734 207115 
1050000 3072049 211314 
1100000 3288502 216454 
1150000 3509106 220603 
1200000 3733965 224860 
1250000 3962171 228206 
1300000 4194158 231987 
1350000 4429331 235173 
1400000 4668560 239229 
1450000 4910513 241953 
1500000 5155601 245088 
1550000 5404652 249051 
1600000 5656512 251859 
1650000 5911306 254794 
1700000 6169686 258380 
1750000 6431383 261697 
1800000 6696172 264789 
1850000 6963206 267034 
1900000 7232974 269768 
1950000 7505412 272438 
2000000 7779763 274351 
2050000 8056579 276816 John McCarthy 
2100000 8336442 279863 17 Mount Street 
2150000 8620053 283611 Mansfield 
2200000 8905641 285588 Notts. 
2250000 9194727 289086 NG19 7AT 
2300000 9486449 291722 United Kingdom 
2350000 9780105 293655 
2400000 10076920 296815 
2450000 10375202 298282 


2500000 10676383 301181 31 


